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Abstract 



We study the "paperclip" model of boundary interaction with the 
topological angle 6 equal to tt. We propose exact expression for 
the disk partition function in terms of solutions of certain ordinary 
differential equation. Large distance asymptotic form of the par- 
tition function which follows from this proposal makes it possible 
to identify the infrared fixed point of the paperclip boundary flow 
at 9 = tt. 
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Field theory models where interactions emerge as a consequence of geo- 
metric constrains imposed on the fields are much interest. They are esthet- 
ically attractive, and they often have rich content and important applica- 
tions. Typical models of this class are nonlinear sigma models where the 
constraints are imposed in the bulk of the space-time. Lately the other type 
of models, where the constraints are imposed only on a boundary, attracts 
much attention. Two-dimensional models of this type emerge naturally in 
string theories and in some condensed matter problems. In string theories 
they provide the world-sheet description of "branes", while in condensed 
matter theory they describe either quantum impurities or "quantum dots" . 

In this paper we study one model of this type. It is the so-called paperclip 
model introduced in [1]. This two-dimensional model of quantum field theory 
involves two-component Bose field X(z, z) = (X(z, z), Y(z, z)) living on the 
disk of radius R. In the bulk, i.e. at \z\ < R, the field X(z,z) is a free 
massless field, as described by the bulk action 

Auik[x] = - / d 2 z d z x-d 2 x, (i) 

* J\z\<R 

while the boundary values Xg of this field, X# = X|| 2 | =fi , are subjected to 
a nonlinear constraint 

r cosh (^f) - cos (^)=0, \ Yb \<^V^+2. (2) 

Here n and r are real and positive parameters (despite the notation, n is 
not necessarily integer). The renormalization does not affect the parameter 
n, which is thus a scale-independent constant, while r "flows" under the RG 
transformations; up to two loops, the flow is described by the equation 

K = (n+1) (l-r 2 )r n , (3) 

where k = %■ is inversely proportional to the RG energy scale E. Here E* 
is the integration constant of the RG equation, which sets up the "physical 
scale" in the model. As in [1], we will always choose E equal to Br 1 , the 
inverse radius of the disk, so that 

k = E*R. (4) 

The equation @ defines a closed curve in the (Xb,Yb) plane, which at 
sufficiently small r has a paperclip shape (see Fig.^), hence the name of the 
model. 
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Figure 1: The paperclip shape for r = 2 x 10 3 . 

When n goes to infinity and simultaneously r goes to 1 as r — 1 — 2(n+i)g 
with finite g, the paperclip curve becomes a circle 

Xi + F| = l, (5) 
and the flow equation (J3J) reduces to 

k = g' 1 e"^ . (6) 

The boundary constraint © defines the "circular brane" model [2]. 

The "circular brane" model has important application in condensed mat- 
ter physics. It is equivalent to the model of dissipative quantum mechanics 
known as Ambegaokar-Eckern-Schon model [3]. The latter is used to de- 
scribe the low-energy sector of the model of a "weakly blockaded" quantum 
dot introduced in [4,5] (see also the review article [6]). That model describes 
an almost open dot which has a large number of energy levels and a large 
charging energy. The dot is connected to the bulk via n degenerate chan- 
nels, and the Ambegaokar-Eckern-Schon model applies when the number of 
channels is very large. In this application the circumference 2ttR of the disk 
in has obvious interpretation as the inverse temperature 1/T, while the 
large charging energy provides explicit ultraviolet cut-off for the model, and 
the topological angle of the circular brane (introduced later in the text) is 
related to the gate voltage. Although the Ambegaokar-Eckern-Schon model 
appears only in the n — > oo limit of the model (J2J) when the paperclip be- 
comes the circle (|5*|). we have reasons to believe that the paperclip model 
with finite integer n has certain relation to the quantum dot model with 
n open degenerate channels. We intend to explore this possible relation in 
separate work. 

The primary object of interest is the boundary state \B), in particu- 
lar its overlap with the Fock vacuum | P ) characterized by the zero-mode 
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momentum P = (P, Q) of the free field X = (X, Y). The overlap 

(P\B) = Z(P\k) (7) 
can be expressed through the disk one-point function 

i? i/3-PV2 z(P\k) = ( e ip ' x (0, 0) > d . sk = J DX e ip - x (°'°) e~^[x] , (8) 

where the functional integration is over all fields X(z, z) obeying the bound- 
ary constraint ©• 

Let us repeat again that Eq. © , as well as the paperclip equation (J2J) , was 
obtained perturbatively, in the two-loop approximation. Therefore it pro- 
vides useful description of the boundary condition only in the weak coupling 
regime, where the curvature of the paperclip curve © is small everywhere; 
this requires n to be large, n> 1, and r to be sufficiently small (so that 
r n <C 1); according to ©, the last condition is fulfilled at sufficiently small 
R, therefore Eqs. (|2j), © provide ultraviolet (UV) description of the bound- 
ary condition. At large distances (large R) and at n ~ 1 the higher loops 1 
and non-perturbative corrections are important. Generally, one expects that 
at R — > oo the paperclip boundary theory "flows" to some infrared (IR) fixed 
point, and so the R = oo limit of the boundary state | B ) is described in 
terms of conformal boundary theory associated with the IR fixed point. 

In regard to the non-perturbative effects, it is important to realize that 
under general definition the paperclip model may involve an additional pa- 
rameter, the topological angle 8. Since topologically the paperclip curve © 
is a circle, the configuration space for the field X(z, z) consists of topological 
sectors, each characterized by integer w which is the number of times the 
boundary value X# winds around the paperclip curve when one goes around 
the disk boundary \z\ = R. This allows one to add the weight factors e lwS 
to all contributions to the functional integral coming from the sectors with 
the winding number w. Thus, in general 

oo 

Z(P|«)= e iwe Z^(P|K), (9) 

w=— oo 

where receives contributions from the topological sector w only. Of 
course, the contributions of the instanton sectors w/0 are invisible in the 

1 The higher loop corrections to @, © are scheme dependent. There are reasons to 
believe that Eqs. (0,© are perturbatively exact, i.e. a scheme exists in which these 
equations are exact to all orders in the loop expansion, see [1]. 
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perturbation theory, and therefore the UV limit R —* of the theory is 
insensitive to the topological angle (this statement is not limited to the case 
n>l; the instanton contributions are suppressed by the powers of k at any 
re, see [1]). However, the IR behavior of the theory can depend on 9 in a 
significant way. 

The paperclip model is in a close analogy with the so called "sausage" 
sigma model proposed and studied in [7] (perhaps better known is the sym- 
metric limit of the sausage model, that 0(3) nonlinear sigma model; its 
counterpart is the circular brane model ©)• The IR physics of the sausage 
(and its 0(3) limit) sigma model strongly depends on its topological angle 
9 [8,9]. The sausage sigma model is believed to be integrable at two val- 
ues of the topological angle, 9 = and 9 = tt. The arguments were given 
in [7], where exact solutions of the sausage model at these two values of 
9 were proposed. According to this proposal, the sausage model at 9 = 
is massive and its solution is described in terms of factorizable S"-matrix 
of three massive particles. On the contrary, the sausage model at 9 = ir 
"flows" in the IR limit to a critical fixed point, which is the c = 1 CFT of a 
compactified free boson; the "flow" is described in terms of certain massless 
Thermodynamic Bethe Ansatz (TBA) system [7]. Solutions for the 0(3) 
sigma model at 9 = and 9 = tt were previously proposed in [10] and [11]. 
By the analogy, one may expect that the paperclip model is also integrable 
in two cases 9 = and 9 = tt. 

There is more then just this analogy to suspect that the paperclip model 
might be integrable at least at some values of 9. In [1] a set of commuting 
local integrals of motion { I 2 i— l }fZi °f the free Bose theory Q was displayed 
which has all symmetries of the paperclip (|2j), and moreover it was found 
that given the symmetries the set is essentially unique. The idea that this 
series has something to do with the paperclip model can be supported by 
the analysis of the UV properties of the model. In the UV limit R — > the 
parameter r in (J2J) becomes small, and the paperclip grows long in the X 
direction; in this limit the paperclip can be regarded as the composition of 
left and right "hairpins" (see Fig.[2J). 




Figure 2: The paperclip formed by a junction of two hairpins. 
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By the "hairpins" we understand the curves 



|exp(±^)=cos(^=), (10) 

where the sign plus (minus) has to be taken for the left (right) hairpin. If any 
one of the curves (|10|) is taken as the boundary constraint, the associated 
left and right "hairpin models" are conformally invariant, and moreover each 
has extended conformal symmetry with respect to certain VF-algebra [1]. 
Although the ^-algebras of the right and left hairpin models are isomorphic, 
the generators of W^ c ' and W^' are realized by different operators in the 
space of states of the free boson theory However, the sets of operators 
iy( c ) and have nontrivial intersection, which is exactly the set of 

commuting integrals {l2z_i} [1]. For these reasons we refer to the set {I2Z-1} 
as the "paperclip series" of local IM. 

In [1] an exact expression for the amplitude Zg = o(P\ k) was proposed 
in terms of solutions of certain linear differential equation. The proposal 
was inspired by remarkable observation of Dorey and Tateo [12] who found 
that in somewhat simpler integrable model of boundary interaction (the 
minimal CFT with non-conformal boundary perturbation) the overlap am- 
plitude analogous to (Q) is related to certain monodromy coefficients of the 
Schroedinger equation with the potential \x\ 2a . Since [12], this finding was 
confirmed and extended to several other integrable models of boundary in- 
teraction [13-15]. Although it is fair to say that true roots of this relation re- 
main mysterious to us, the examples suggest that the relation may be rather 
general. Given a model of boundary interaction suspect of being integrable, 
it is worth trying to identify associated ordinary differential equation. The 
proposal of [1] was made according to this strategy. Let us briefly summarize 
it here. 

The main ingredient is the ordinary differential equation 



d2 nP 2 eX _ (n+2)Q 2 -i e x 2 , s 

dx 2 4 1 + &x 4 (1 + e ^2 + < I 



= 0, (11) 



where P and Q are the X and Y components of the zero-mode momentum 
P, and k is the same as in @; according to the definition @ we assume for 
the moment that k is real and positive. 

Eq. (jllj) has a form of a stationary zero energy Schroedinger equation 
with specific potential V(x) given by the last three terms in (|llj) . The 
potential V(x) is positive and grows fast at large positive x, therefore (jTTj) 
has a solution H(x) decaying at x — * +00; this condition specifies E(x) 
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uniquely up to normalization 2 . To fix the normalization, we assume that 

S{x) - exp | - (f + k) x - K J* ^ ((1 + - 1) 1 (12) 

as x — > +oo. On the other hand, V{x) approaches positive constant k 2 
at large negative x. Hence Eq. (jllj) has a solution which decays at large 
negative x; we denote this solution ^ + (x). The condition, 

* + (x) — ► — r as x — > oo , (13) 

+ v ; r(l + 2«) V ; 

specifies the solution \I'+(x) uniquely, including its normalization. Then 

_ /2k\ 2k 

Z 9=0 (P\k) =g 2 D V^ (— ) ^[2,*+]. (14) 

Here and below gr> = 2~ 1//4 is the ^-factor [16] of the Dirichlet boundary, 
and W[F, G] denotes the Wronskian F(x)G'(x) - F'{x)G(x). Eq. (|T1|) is the 
proposal of [1]. 

It was shown in [1] that (|14[) exhibits the UV (i.e k — > 0) behavior com- 
pletely consistent with what one expects from the "hairpin decomposition" 
of the paperclip at R — ► 0. At small k the expression Q14|) can be written as 
(here we write explicitly the components of P = (P, Q)) 

Z e=0 ( P, Q | «) = B{P, Q) P e=0 (P, Q | k) + B(-P, Q) F e=0 (-P, Q\k), (15) 

where 

B(P,Q)=& (16) 

v^r(-i^p)r(i-iP/^) 

rQ- 2-i^f)r(i + v / ^+2 § -iv^f) ' 

and F$ = o(P, Q \ k) (apart from the factor k 2k ) admits asymptotic expansion 

2 

in a double series in powers of k and k> (here and below we use the symbol 
~ to indicate relations which hold in the sense of asymptotic series), 

oo 

Fe=o(P, Q\k) ~ k 2k Ai( P > ^ + " . ( 17 ) 

i,j=Q 



2 Here we slightly change notations for the solutions of as compared to [1] 
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with /o ; o = 1- From (|14|). the powers appear from the k — > expansion 
of the solution while the integer powers of k come from the expansion 

of ^ + (x). 

The expansion (|15|) is in good agreement with expected form of (jSJ) in the 
domain «< 1. Recall that according to (@J) and (jSJ this domain corresponds 
to the UV limit of the paperclip model, where the paperclip © can be 
regarded as the composition of two hairpins, as in Fig.|2j Roughly speaking, 
the two terms in (|15|) correspond to contributions from the right and the left 
hairpins, respectively. More precisely, at small n and QmP ^ the saddle 
points of the functional integral Q correspond to the field configurations 
where the boundary values (Xb, Yb) are mostly concentrated near the left 
or the right end of the paperclip ©. The factors B(P,Q) and B(—P,Q) 
are precisely the boundary overlap amplitudes of the right and left hairpin 
models respectively. The powers of k« in (j!7[) are associated with 

the perturbative corrections, due to "small" fluctuations around the saddle 
points, which feel only the small deviations of the shape of the paperclip 
© from the respective hairpin, while the integer powers of k are associated 
with the contributions of the topological sectors in Q (see Ref. [1] for 

more detailed discussion, which also includes the explanation of the factor 
k 2k in (|14|) and (|17|) as the effect of small instantons). 

On the other hand, when k — > oo (and P, Q are fixed) the potential term 
in (jllj) becomes large at all real x, and therefore the Wronskian in (|14|) can 
be found by straightforward application of WKB technique. This results in 
the asymptotic k — > oo series 

Z e=0 (P \ k) c g 2 D eX p{-f;%^H, (18) 

i=i K J 

where I21-1O?) = ^2l~i{P-,Q) are certain polynomials in P 2 , and Q 2 of 
the degree I, with coefficients which depend only on n. It turns out that 
these polynomials are exactly the vacuum eigenvalues of the local IM of the 
"paperclip series" {^-l}, 

|P> = R 1 ' 21 hi-i{P,Q) |P> (19) 

(see [1] for details). This form of IR expansion suggests that at 9 = the 
IR fixed point is just the Dirichlet boundary condition for the free field 
X = (X,Y): 

(X B ,Y B ) = (0,0). (20) 
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According to Q18|) the RG flow approaches the Dirichlet fixed point along 
irrelevant direction which is a combination of densities of the paperclip IM 
{hl-i}- 

In this paper we extend the proposal of [1] to the case 8 = tt. The 
starting point is the same differential equation but now instead of ^ + 
in (|14|) . we take another solution which grows as e~ KX at large negative x. 
Of course, this condition alone does not define the solution uniquely since, 
besides overall normalization, one can always add any amount of *$>+(x). 
Usually it is difficult to make unambiguous definition of a growing solution, 
but in our case the following property of helps. Let us consider x 
as a complex variable. The potential V{x) is analytic function of x with 
the branching-point singularities at all points where e x turns to —1. Let us 
make brunch cuts from each of the points x = iri (2iV + l), N = 0, ±1, ±2, . . . 
to +oo parallel to the real axis, and choose the branch of V{x) on which 
(l + e x ) n is real and positive on the real axis of the x-plane. Now, restricting 
attention to the domain ?R.ex < one finds that the potential V(x) has the 
periodicity property 

V(x + 2vri) = V{x) (Ke x < 0) . (21) 

Consequently, the equation (jllj) has two Bloch-wave solutions (2k ^ Z): 

V±(x + 2TTi) = e ±2wiK y±(x) Oftex<0), (22) 

where the Bloch factors are found by taking the limit $tex —* — oo. At this 
point we assume that 2k is not an integer, so that the conditions (|22j) specify 
the two independent solutions ^?±(x) uniquely, up to their normalizations. 
Of course, the solution ^+(x) defined this way decays as e KX at 3?e x — > — oo, 
and the asymptotic condition also fixes its normalization. The solution 
^f-(x) grows at large negative Kex, and its normalization can be fixed by 
specifying the leading asymptotic in this domain. Thus we define (x) by 
the conditions 



^„(x + 2vri) =e~ 2 ^ iK ^_(x) (Kex<0), (23) 

q KX 

^-(x) —>■ — r as $iex^— oo . 



It is possible to show that both ^+(x) and ^_(x) defined by and (|23|) 
are entire functions of k, and 

*_(«|a;) = a?), (24) 
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where we temporarily exhibited the dependence of *$>± of the parameter k. 
From the definitions and Q23JI we have 

W\V-,*+ = — i -. (25) 

Our proposal for Zg =n is 

Z e=7T (P\ K ) =g 2 D V^ (^)" 2 V[H,*_] . (26) 

The k — ► form of the Wronskian in Q26|) can be derived through the 
perturbative evaluation of the solutions and ^_(x), as it was done 

in [1] for the Wronskian in (|14|) . This leads to the expansion similar to 

CESJ), 

2*=*(P,Q|«) =B(P,Q) F e=7T (P,Q\ K )+B(-P,Q) F e=n (-P,Q\ K ), (27) 
where P(-P, Q) is the same as in (|15j) . while 

oo 

F e=7r (P, Q|«) ~ kT 2k ^ /^-(P, Q) (-«)* K v , (28) 

i,j=Q 

with exactly the same coefficients fij(P,Q) as in (|17|). This form follows 
from the fact that integer powers of k in (|28j) come from perturbative expan- 
sion of the solution ^^(x); in view of (|24|) they are related to corresponding 
integer powers in (|17j) by the change of the sign, k — > — k. At the same time, 

2 

the powers of k« appear as the result of expansion of H(x), and hence they 
remain unchanged in (|28|) as compared to (|17|) . Since the powers of fi> are 
interpreted as the perturbative contributions to the functional integral (jSJ), 
while the integer powers of k are due to the instanton contributions, the 
form (|27|) . (|28|) is exactly what one expects to have from the definition Q 
at 9 = ir. This property of the UV expansion was the main motivation of 
our proposal (|2^)l . 

Accepting (|26j) . we can address the problem of the IR behavior of the 
paperclip model at 9 = ir. One has to find the k — > oo asymptotic of 
the Wronskian in (|26|). Unlike the case = 0, Eq. (|14|) . this turns out to 
be rather subtle problem. While at k — > oo the WKB approximation still 
formally applies to the problem is analogous to the problem of finding 
the "over-the-barrier" reflection amplitude in quantum mechanics which in 
the semiclassical approximation requires identifying appropriate "turning 
points" (the zeroes of the potential V(x)) in the complex x-plane [17]. In 
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our case, at large k all the complex turning points approach the singular 
points e x = — 1 where the semiclassical approximation for breaks down. 
Finding correct k — > oo expansion of (|26|) requires analysis of the solutions 
of in the vicinities of the singular points e x = —1. The calculations 
are rather involved, and will be present them elsewhere. The result is the 
asymptotic k — > oo expansion 

Z e=n (P\ K )^g 2 D ^ T(P\k) exp{ V ^zlg) 1 (29) 

L i=i K J 

where a is a constant 3 , i2Z-i(P) are the same eigenvalues of the "paperclip 

IM" I2Z-1 as in (|18j) . and T(P | k) is the asymptotic series in inverse powers 

2 

Of K n + 2 , 

00 

T(P|k)~2cos(^=)+J>(^<3)«~^- (30) 

Z=l 

Here again P and Q are t ne an d V components of the zero-mode mo- 
mentum, i.e. P = (P, Q), and the coefficients ti(P,Q) are in principle com- 
putable through perturbative solution of (|11|) in the vicinity of the singular 
points e x = — 1. The first two coefficients can be evaluated in closed form. 
We have 



re + 2 (n + 2)Q 2 -nP 2 -1 
n + 4 2((n + 2)Q 2 - 1) 



27T a 



w_y la p/" _ _y la 



while the expression for t2(P,Q) is somewhat cumbersome, and we present 
it in Appendix. 

The most important difference of (|29|) from Q18[) is in the factor (|30|) 

2 

which involves the powers of ft «+ 2 ; its presence indicates that in this case 
the IR fixed point differs from the trivial Dirichlet boundary constraint 1)20(1 . 
The series ()30|) starts with the term 2 cos ("^7=^5) which suggests that in the 
limit k — > 00 the boundary values of X are constrained to two points, 

(X B ,Y B ) = (0,± 7 ^). (32) 



3 a = 2 7B + 2^(l + f) with i/>(z) = & lo g r(«), 7e = 
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Eq. (|32j) characterizes the IR fixed point of the paperclip boundary flow at 
9 = it. At large but finite n the boundary value (Xb, Yb) is allowed to jump 
between the two points Q32JI . the possibility of such jumps being responsible 
for the higher-order terms in the series (J30J). In other words, the RG flow 
arrives at the IR fixed point along certain irrelevant boundary fields 
which generate such jumps. Obvious candidates are the boundary vertex 
operators 

V ± {r) = [d a X] B ^^(t), (33) 

where r is natural coordinate along the boundary (we write z/R = e~fe~, 
so that the boundary \z\ = R is at a = 0), [d^X]s = d a X(a, t)| ct= o is the 
normal derivative of the field X at the boundary, and Yb(t) is the boundary 
value of Y, the T-dual of Y 4 . The exponentials exp (± ^^ 2 ) (t) create 
jumps in the boundary value Yb at the point r, i.e. the boundary values Ye 

/n+2' 



to the right and to the left of each vertex V± (r) differ by the amount db- 

Note that the dimensions of the fields (|3*3*|) equal to 1 + > which is in exact 

agreement with the fractional power k,~™+z appearing in the expansion (|3U|) . 
According to these arguments, it seems likely that the full asymptotic series 
Q3UJI can be generated by expanding the following expression 

T(P|k)= (34) 

r j^-Q r f 2nR 1 1 

(PlltJevW TexpjA y dr (cr_y + (r) + ct+F_(t)) |J | 5 D > . 

Here | -Bd ) is the boundary state of the free theory © with the Dirichlet 
boundary condition (|2U|) which is the superposition | Bjj ) = f d 2 P | Ip ) 
of the Dirichlet Ishibashi states \Ip). In what follows we assume that the 
exponentials in (|33|) are normalized in such a way that 

(P\V+(t)V-(t')\B d ) -» 2 Ir-rp 2- ^ as r-r'^0. (35) 

The r-ordered exponential is understood in terms of its expansion in powers 
of the parameter A (which will be related to the parameter , see Eq. (|3*7|) 
below), 

r(P|K)=2cos(^) + f j 2 J R &r j\r'x (36) 
"Tifs F + ( T )y_(r') + evlfj y_(r)y + (T')l | > + • • • • 



4 The T-dual of the free massless field is defined as usual, through the relations d T Y — 
idvY, and d a Y = -id T Y . 
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The relations a\ = guarantee that the boundary values of Y are limited 
to two points indicated in (|3*2*)) . With this understood, two eigenstates of a z 
with the eigenvalues ±1 appear to be in correspondence with two possible 
boundary values Yb = ±-^=. At the same time, the a-matrices in (|34"|) 
can be regarded as operators representing an additional boundary degree 
of freedom; then the expression in the exponential in (j34j) is interpreted 
as the boundary action describing perturbation of the Dirichlet boundary 
condition. Anyhow, the g-factor at the IR fixed point (i.e. (|29|) evaluated 
at k = oo and P = 0) is 2g 2 - ) , not g 2 ^ as it was in the case of the 6 = IR 
fixed point (|20|) . 

As was mentioned above, the dimensions of the vertex operators (j33|) are 
greater than one. As the result, the integrals appearing in the expansion 
(I36|) generally diverge when r-separations between the insertions V±(t) be- 
come small. One has to specify some regularization to make the expression 
(jSU meaningful. We assume here the "analytic regularization" common in 
conformal perturbation theory. The prescription is to consider n as an arbi- 
trary complex number and to evaluate the integrals in the domain of n where 
they converge (specifically, at IRere < —4), and then to continue in n to the 
real positive values. It turns out that this continuation yields unique finite 
results for all integrals involved in the expansion (HU). Thus, the A 2 term 
explicitly written in (|36j) is readily evaluated in analytic form. Remarkably, 
its dependence of the momenta (P, Q) appears exactly the same as in (|31j) . 
and upon identification 



it coincides with the first (/ = 1) term ti(P, Q)k~™+2 of the expansion 
(jSOJ). We believe that under the identification ()37j) . and with the analytic 
regularization described above, the expression ()34|) reproduces all terms of 
the asymptotic series (|3*U)l . i.e. it plays the role of the "effective IR theory" 
describing the approach of the 9 = ir paperclip flow to its IR fixed point 
()32|) . Note that the effective theory (|3*1)) is expected to capture only the 
"perturbative" part (|3*U|) of the full IR expansion (|29|) . and it does not seem 
to have any control over the "nonperturbative" terms due to the local IM 



The relation (|37j) is singular at n = oo. Correspondingly, obtaining 
correct n — > oo limit from the effective theory ()34|) is not exactly straight- 
forward. At n > 1 alternative form of the IR effective theory is more 




(37) 



2 



in (JH. 
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convenient. It is written as 



r(P|«)= (38) 

(P | TrTexp { J dr a x [d a X] B (r) + ^ a y [d a Y] B (T))j\ B D ) , 

where again and are the normal derivatives of the fields X 

and Y at the boundary, and a x , a y are coupling constants. Formally, (|38jl 
can be brought to the form (|34|) by field-dependent gauge transformation 

ffa - e"^ (r) a a e 1 ^ (r) , (39) 

with obvious relabeling of the c-matrices to match the notations in (|34|). The 
subtlety in this transformation is in the renormalization of the parameters. 
Perturbation theory in the couplings a x , a y in (|38jl has logarithmic UV 
divergences which lead to renormalization of these parameters. In view of 
the X <-> y symmetry the corresponding RG flow equations can be written 
as 

-E^ = P(g x ,g y ), -E^ = (3(g y ,g x ), (40) 

where E is the RG energy scale, and g x = a x , g y = a y . The model (|3*S|) 
was studied in [19], where the leading (one-loop) term of the beta-function 
is presented, 

f3(g x ,g y ) = -2 g x g y + 0(g 3 ) . (41) 

The higher loop terms depend on the renormalization scheme. One can note 
that l)34|l is easily evaluated in the limit when g x — > with g y kept finite 
(or when g y — > with finite <7 X )- It follows from this solution that a class of 
"natural" schemes exists in which 

(3(g,0)=f3(0,g) = 0, ±p(g,h)\ g=0 = -2h. (42) 

In any such scheme 

P(g„ Sy) = "2 g x g y + b g\g y + 0(g 4 ) , (43) 

where the constant b is scheme independent. Explicit two-loop calculation 
yields 6 = 4. The terms ~ g 4 and higher still depend on the scheme. It 
looks likely that any RG flow equation of the form (|4U|) with the /3-function 



13 



satisfying (|42|) and having the expansion (|43|) with b = 4 can be brought, 
order by order in g, to the following convenient form 

-E^ = -E^ = (44) 

d£ d£ 1 + 5x + 5y 

where g x = <7x/(l — <?x)> <?y = 5y/(l — 9y)> by appropriate redefinition of the 
coupling constants <? x , g y (we have explicitly checked this statement up to 
the order g 6 ). The RG flow ()44l) conserves the difference e = g y — g x . In the 
coordinates g x ,g y the RG trajectories are straight lines which end at the IR 
fixed points g x = if e > 0, or g y = if e < 0. Let us assume here e > 0. 
The RG invariant e can be related to the parameter n in (|34j) . 

e = 9y ~ 9x = — r • (45) 
n + 1 

Indeed, at the IR fixed point g x = while g y = the last quantity must 
coincide with the dimension of the coupling constant A in ()34|) . Integrating 
(jSJ) one finds 

9 * = — ~p' g ^ = —p^ (46) 
where p decays towards IR according to the equation 

€ph - E (47) 



1 - P £o 

where Eq is the integration constant; one can set Eq = C E* with an arbi- 
trary constant C. The equations (|1d]> . (|4*7|) define a family of ^-dependent 
renormalized coupling constants g x , g y . It is clear that (|38|) admits pertur- 
bative expansion in these coupling constants, and if one chooses E = R^ 1 
the coefficients of this series are undependent of R. 

Let us consider the case n = oo which is of special interest because it 
corresponds to the circular brane model [2] at 9 = tt, and hence is equivalent 
to the Ambegaokar-Eckern-Schon model [3] of quantum dot at the degen- 
eracy point. Obviously, in this limit the form (|3U|) of the IR asymptotic 
expansion of the boundary amplitude is no longer meaningful. Instead, in 
this case the factor T{ P | k ) in ()29l) rather expands in powers the "running 
coupling constant" g related to k by the equation 

9 = -zr ■ ( 48 ) 

e c ° k 
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The same relation is obtained by integrating the RG equation (|44j) with 
9x = 9y = 9 ■ It is possible to show that at n = oo the factor T( P | k ) in 
(|29j) expands as 

T( P | k) = 2 - (vrP) 2 g + ( vr 2 - 2 (^P) 2 + i (^P) 4 ) g 2 + 0(g 3 ) , (49) 
where g = g{n) is determined through Q48[) with 

co = 3 + 7 i? , (50) 

which is in perfect agreement with the results of renormalized perturbation 
theory in (|38|) . 

The differential equation can be integrated numerically, thus pro- 
viding high precision numerical data for the disk partition functions (|14|) 
and (|26j) . To handle the circular-brane case n = oo an appropriate limit of 
the differential equation (jllj) has to be taken. Making shift of the variable 
x = y — logn and then sending n to infinity, one obtains the differential 
equation (124) of Ref. [1], while the asymptotic condition for the solution 
H is also suitably modified, see Eq. (126) of Ref. [1]. 

In Table ^ we compare data for the boundary entropy [16, 20] of the 
circular brane (Ambegaokar-Eckern-Schon) model 

S 6 {k) = (i-k^) log^(0|«) , (51) 

obtained by numerical integration of the differential equation corresponding 
to n = oo with first few terms of the IR expansions 

AS^o(k) = J---J-3+ J 1 5 +0(k- 7 ), 
6 k 180 k 105000 k° 

AS^(k) = log(2) + ^5 2 + 2vr 2 5 3 -^^5 4 + 0(5 5 )- (52) 

1 5" z+0(K ') , 

6k 180 k 3 105000 k 5 v ; ' 

where ASq = Sq — log(g| ) ). The boundary entropy for 9 = 0, ir as a function 
of the dimensionless inverse temperature k = A, is plotted in Fig.Ol 

Omitting details, we would like to mention here that using methods de- 
veloped in [21] (see also [13,23]) one can reduce the calculation of the Wron- 
skians (|14|) . (|26jl at special values of the parameters to solution of certain 
system of nonlinear integral equations, the "TBA system". For example, for 
integer n and 

P = 0, Q = ±_L=, (53) 
V n + 2 
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0.04749 


0.8291 
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0.84827 


4.0 


0.0416 


0.04158 


0.8231 


0.83989 


4.5 


0.0370 


0.03698 


0.8179 


0.83298 


5.0 


0.0333 


0.03329 


0.8136 


0.82716 


5.5 


0.0303 


0.03027 


0.8097 


0.82216 


6.0 


0.0278 


0.02775 


0.8064 


0.81780 


6.5 


0.0256 


0.02562 


0.8034 


0.81396 


7.0 


0.0238 


0.02379 


0.8006 


0.81053 


7.5 


0.0222 


0.02221 


0.7980 


0.80745 


8.0 


0.0208 


0.02082 


0.7959 


0.80466 


8.5 


0.0196 


0.01960 


0.7938 


0.80212 



Table 1: Comparison of the boundary entropy ASe=o j7T = £0=0,77 — log(g|>) 
for the circular brane as computed from the differential equation with few 
terms of the IR asymptotic expansions AS^ 0n given in Eqs. (|52j) . 

the amplitude Z(P | k) (J7| is expressed through solution of the D n -type 
TBA system, associated with the diagram in Fig.0J It is remarkable that 
this system coincides exactly with the system of TBA equations for the 
Toulose limit of the n-channel Kondo model [22] , while the entropy g( Kondo ^ 
of the impurity spin j = is given by 

^(Kondo) = K ^ ) tog ( g~ 4 P* | Z 0=7r (P* | «) ) , (54) 

where P* = (0 , ) and is identified with the Kondo temperature. 

Notice that at this choice of P, Q the coefficient t\ in the low-temperature 
expansion of the partition function (|30|) vanishes and the expansion starts 
with t2 (see TableElin Appendix). Also, this TBA system (see Fig.EJ) differs 
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ASf 




K 



Figure 3: The boundary entropy ASq = S$ — log(g£>) (|51j) for the circular 
brane model, as a function of the dimensionless inverse temperature k for 
9 = 0, 7r, derived from (JTIJ) and (|2oT) . 



only in the structure of the "source terms" from the TBA system associated 
with the Hn \ models [24], which are certain integrable perturbations 
of the minimal Z„-parafermionic CFT. 



1 




n-1 



o 



Figure 4: Incidence diagram for the TBA system in the case of the paperclip 



brane with an integer n > 3 and P = 0, Q 
indicated near the corresponding node. 



/n+2 ' 



The source term 



At n = oo, the boundary entropy of the circular brane model with 6 
0, 7T is given by 



AS e=k7T (K) 



dp 
2^ 



,7-/3 



log [1 + 



-e*Gf3)l 



cosh (f3 — 7) 



(k = 0, 1) 



(55) 



where 7 = log(47TK) and k = 0, 1, 2, . . . solve the infinite chain of integral 
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equations: 



, , f°° d(3 log [(1 + e-^wC)) (1 + e - £ *-i^))] 

= e k (a) + / — — h 

J_ 00 2tt cosn(p — a) 

r^ MH^I (t . 2 ,3,4...c o) , (50) 
7_ 00 2vr cosh(/3 - a) 

f°° dB log [1 + e~ £2 ^] 

e a S kfi = e k (a)+ ^ -^J- -i (fc = 0, 1) . 

7^ 2vr cosh(/3 - a) 



To sunmmarize, we have proposed exact partition function of a nontrivial 
model of boundary interaction - the paperclip model. There are parallels 
between this model on one side, and the bulk sigma-models such as the 
0(3) nonlinear sigma model and the "sausage" model [7] on the other. In 
particular, according to our proposal, the addition of the topological term 
with 9 = ir in the paperclip model leads to a non-trivial IR fixed point. In the 
limit when the parameter n goes to infinity the paperclip model reduces to 
the "circular brane", which is equivalent to the Ambegaokar-Eckern-Schon 
model of a quantum dot. At finite integer n our proposal suggests intriguing 
formal relation to the n-channel Kondo model in the Toulose limit. 



Acknowledgments 

The authors grateful to Lev B. Ioffe, Mikhail V. Feigel'man for valuable 
discussions and Dmitri M. Belov for help with graphical software. SLL and 
ABZ are grateful to Vladimir V. Bazhanov, Vladimir A. Fateev and Alexei 
B. Zamolodchikov for sharing their insights. AMT is grateful to Igor L. 
Aleiner for discussions on quantum dots. 

The research of SLL and ABZ is supported in part by DOE grant #DE- 
FG02-96 ER 40949. SLL also acknowledges the support from the Institute 
for Strongly Correlated and Complex System at BNL, where most of the re- 
sults of this paper have been obtained. AMT acknowledges the support from 
US DOE under contract number #DE-AC02-98 CH 10886. ABZ gratefully 
acknowledges hospitality and generous support from Foundation de l'Ecole 
Normale Superieure. 



18 



Appendix 

Here we present an integral representation for the second coefficient ap- 
peared in the expansion (|3U|) . which was obtained by an analytic technique 
developed in Appendix A of Ref. [18]. Bellow we use the following notations 



p 



P , Q 

n — , q = yn + 2 — 



(57) 



n + 2 1 1 v 2 
The coefficient t 2 can be written as 

t 2 = -4 sin 2 «) cos (2tt£(<? + 1)) a\ - 4 sin(27r£) sin (2vr(g + a 2 • (58) 

Here 



oi = (20 



_ 2g r(Qr(i-Q 

2^ 



1 - 2£ 



+ 



r(i + (2 g + i)Q 
r((2?-i)£) 



,(59) 



2(1 + 2^ 2(g2_i; 
and a2 in (|58|1 admits the following representation for $te q > and < £ < 



«2 



r(i + 2g( g + i)) / a , R , r \ 
r(2^-i)) l^T + i? + C )' 



(60) 



where 
A 



B 

C 



(20 



0Fr(i-2Qr 2 (i + 3Q 
4£r2(i + £)r(i + 2e) 
r(i + 2g)r(|-2g) 

8^(1 + 40 ' 



-44 



-4£ 



1 — 2^ 2p 2 
2(1 + 20 + ~3" 



(61) 



8r2(i-or2(i + o 

dr F(r) sinh(27rr) 



2vr ( r + i^)( r -i^)(r + i0 



1 — 2$ 



2 ^2 



2^p 



2(1 + 20 4r 2 + e 2 



with 



F(t) = T(l - 2$ + 2ir)r(l - 2$ - 2ir)r 2 (l + $ + 2ir)r 2 (l + £ - 2ir) . 



19 



Notice that for £ = -r (n = 2) one has to use the L'opital's rule to evaluate ai 
from Eq. (|60|) . Using (|58(1 it is straightforward to study the large n-behavior 
(£-0) of t 2 : 



27T 2 f 



-44 



7T 



1 



2 4 
3? 



2<? 2 - 2p 2 - (3 

-(£-!)p 2 - 



4g 2 + 4p 2 



° 2 2 ' 

3 Q P 



£ + (62) 



The above representation is very useful in numerical evaluation of this 
coefficient because of the fast convergency of the integral in (|61|) . In Table 
2 we present the numerical values of the coefficient at some integer n and 
P, Q given by Eq. (|53|) . These numbers were also confirmed by numerical 
integration of the TBA system depicted in Fig. 2] 



n 


«2 


1 





2 






48 


3 


0.0546105 


4 


0.0661040 


5 


0.0683646 


6 


0.0658731 


7 


0.0613178 


8 


0.0561029 


9 


0.0509101 


10 


0.0460445 



Table 2: Numerical values of the second coefficient t% in the expansion (|30|) 
for P = 0, Q = 



M+2 " 
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